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Inclusive production of D*^(2010) mesons in deep inelastic scatter- 
ing has been measured using e^p and e~p data obtained with the ZEUS 
detector at HERA using integrated luminosities of 16.7 and 65.2pb~^, re- 
spectively. The decay channel D*^ — > D^tt"*" with K~7r~^ and corre- 
sponding antiparticle decays were used to identify D*^ mesons. The D*^ 
cross sections in e'p and e+p interactions agree with NLO QCD predic- 
tions, although the D*^ cross section in e~p is slightly higher than that in 

1. Introduction 

Charm production in deep inelastic scattering (DIS) at HERA has been 
shown in previous studies to be consistent with purely dynamic Boson-Gluon 
Fusion (BGF) production [1-3]. This agreement has now been tested with a 
larger data sample than the previous ZEUS measurements, and e~p as well 
as e~^p cross sections have been calculated. The charmed mesons were iden- 
tified using the decay D*~^ — > D^n^ with K~7r~^ and corresponding 
antiparticle processes, where tTs refers to a low momentum pion accompa- 
nying the D^. The differential cross sections are measured as functions of 

and Bjorken x, defined as Q"^ = —(^ = {k — k')'^ and x = Q'^/{2P ■ q), 
where k and k' are the four-momenta of the initial and final state lepton, 
and P is the four-momentum of the proton. 



2. Kinematic reconstruction and event selection 

The ZEUS detector is described in detail elsewhere [4,5]. The x and 
variables were reconstructed by the S-method, which uses both the 
scattered lepton and the hadronic system measurements [6]. Standard cuts 
were imposed to select neutral current DIS events [7]. The D*^ mesons 
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were selected in the range 1.80 < M{D^) < 1.92 GeV, 0.143 < AM < 0.148 
GeV, 1.5 < priD*) < 15 GeV, and \r]{D*)\ < 1.5. The number of £>*± 
events determined from a 5 parameter fit 

F(AM) = [Pl/V2^ ■ P3] • exp[0.5 • (AM - P2)^/P3] + PA ■ (AM - m^ f^ 

where PI - P5 are free parameters is 1229ib48 in the e~p data, and 4240ib90 
in the e'^p data. The number of D*^ mesons extracted from empirical wrong 
charge background subtraction within the signal region 143 < AM < 148 
MeV, is 1219 lb 58 and 4239 ± 113 in the e~p and e~^p data respectively. The 
AM distributions are shown in Figure 1, for e~p and e~^p data separately. 

3. Study of systematic effects 

The systematic uncertainties on the measured D*^ cross sections were 
determined by changing the selection cuts or analysis procedure. These un- 
certainties are divided into three groups. 

Event reconstruction and selection: 

• The cuts on y^. yjB, 5, and the vertex position were varied [7]. 

• The cut on the position of the scattered lepton in the RCAL was 

raised. 

• The minimum energy of the scattered lepton was raised. 

• The Electron method or Double Angle method [8] was used to recon- 
struct the kinematics. 

D*^ reconstruction: 

• A higher track quality was required (restriction on the polar angle of 
the track). 

• The transverse momentum requirement of the K and tt candidates 
was varied. 

• The signal region for M{D^) and AM were varied. 
Monte Carlo: 

• The acceptance was calculated using HERWIG [9] instead of RAP- 
GAP [10]. 

The overall systematic uncertainty was determined by adding the above 
uncertainties in quadrature. The normalisation uncertainties due to the 
luminosity measurement error, and those due to the D*"^ and D° branching 
ratios were not included. 
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4. Results 

In the kinematic region 1 < < lOOO GeV^, 0.02 <y < 0.08 and the 
selected D*^ region, the cross sections calculated are 

a{e~p e-D*^X) = 10.20 ± QAd,{stat.)lH{syst.) nb, 

(7(e+p ^ e+D*±X) = 8.94 ± {).24.{stat.)l^{syst.) nb. 

The e^p cross section is consistent with that previously published [1], allow- 
ing for the increase in proton beam energy [7], while the e~p cross section 
is slightly higher. Figure 2 shows the differential cross sections as a func- 
tion of and x compared to the NLO calculation implemented in the 
HVQDIS program [11, 12]. This program is based on the BGF mechanism, 
and uses the Peterson fragmentation function [13], with e = 0.035, to hadro- 
nise the charm quark to a D*^. The mass and renormalisation scales were 
set to ^4m^ + Q"^. The hadronisation fraction /(c D*~^) was set to 
0.235 [14]. The boundaries of the shaded band indicate two extreme values 
of HVQDIS predictions, from changing the charm mass between rric = 1.3 
to 1.6 GeV, and using different sets of structure functions, GRV98IIO [15], 
CTEQ5F3 [16] and a ZEUS NLO fit [17]. The NLO calculations based on 
BGF give a good description of the measured D*^ cross section over the 
full range of and x. For > 20 GeV^, the D*^ cross sections in e~p 
and e^p differ slightly, while conventional charm production mechanisms 
contain no charge dependence on the lepton in these interactions. More 
e'p data is essential to investigate whether this is a statistical fluctuation. 
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Fig. 1. Data (solid dots) for AM = 
{Mktttts ~ Mktt) for e^p data above 
and e^p data below. The background 
from wrong charge combinations in 
shown as the filled histogram. The 
solid line shows the result of the fit 
described in the text; the dashed ver- 
tical line indicates the signal region. 



ZEUS 



ZEUS 



-1- 

• ZEUS (pre!.) 1998-99 ep 

ZEUS (prel.) 1999-00 e"p 

1 I HVQDISm^=1 .3-1.6 GeV 

(GRV98, CTEQ5F3, ZEUS NLO) 



.|Q _ • ZEUS (prel.) 1998-99 e p 
o ZEUS (prel.) 1999-00 e> 
- 1 I HVQDiSm^=1.3-1.6GeV 

(GRV98, CTEQ5F3. ZEUS NLO) 

1 r- 



Fig. 2. Differential D*^ cross sections for e~p and e~^p data as a function of 
on the left and x on the right, compared to the NLO QCD calculation of HVQDIS. 
The inner error bars show the statistical uncertainties, while the outer ones are 
the statistical and systematic errors added in quadrature. The boundaries of the 
shaded band for the HVQDIS prediction correspond to the full uncertainty due to 
the charm mass variation and choice of structure function as described in the text. 
The lower portion of each plot shows the ratio of the e^p to e'^p cross sections. 



